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Abstract 



o 



• The elastic and a-production channels for the ^He+^^^Pb reaction are investigated at energies around the Coulomb barrier (£'iab=14, 
16, 18, 22, and 27 MeV). The effect of the two-neutron transfer channels on the elastic scattering has been studied within the 
Coupled-Reaction-Channels (CRC) method. We find that the explicit inclusion of these channels allows a simultaneous description 
of the elastic data and the inclusive a cross sections at backward angles. Three-body Continuum-Discretized Coupled-Channels 
(CDCC) calculations are found to reproduce the elastic data, but not the transfer/breakup data. The trivially-equivalent local 
polarization potential (TELP) derived from the CRC and CDCC calculations are found to explain the features found in previous 
phenomenological optical model calculations for this system. 

•Keywords: Nuclear reaction ^'^^Pb(^He,^He),^''^Pb(^He,'^He), halo nuclei, Continuum-Discretized Coupled-Channels and 

Coupled-Reaction-Channels calculations . 
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1. Introduction 

Since the early nineties, a considerable amount of experi- 
mental data of reactions induced by the Borromean nucleus 
^He on a variety of targets has been accumulated. In the en- 
ergy regime around the Coulomb barrier, and for medium-mass 
and heavy targets, these data show some common remarkable 
features. First, the elastic scattering angular distribution does 
not follow the expected Fresnel pattern, that characterizes the 
scattering of heavy ions at these energies. Instead, the charac- 
teristic Fresnel peak is partially or completely suppressed and 
the angular distribution shows a smooth decrease as a function 
of the scattering angle. Second, these reactions exhibit a large 
yield of a particles |1)-(5I]. This is clearly a consequence of the 
■weak binding of the ^He nucleus, which can be easily broken 
up by the strong couplings induced by the target. 

Heavy ion reactions between stable nuclei can be understood 
in terms of the strong absorption picture This implies 
the existence of an interaction distance (Rsa) such that impact 
parameters below Rsa are dominated by strong absorbing pro- 
cesses (i.e. deep inelastic collisions, compound nucleus, etc), 
whereas for separations larger than Rsa, the nuclei barely in- 
teract, and elastic scattering dominates. Direct reactions, such 
as nucleon transfer or diffractive breakup, take place only for 
peripheral or grazing collisions. However, the analysis of re- 
cent experiments suggest that this picture may be inadequate 
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to interpret collisions involving weakly bound nuclei. In the 
^HeH-^"^Pb case, optical model (OM) calculations performed 
at Coulomb-barrier energies show that, in order to re- 

produce the elastic data, one requires optical potentials with 
a very diffuse imaginary tail. If the interacting potential is 
parametrized using a standard complex Woods-Saxon form, the 
diffuseness parameter of the imaginary component required to 
reproduce the elastic data turns out to be around a,- x 2 fm, a 
value that is considerably larger than the diffuseness derived 
from the matter distribution 0.6 fm). This phenomenon, 
recently referred to as long range absorption effect, suggests 
the presence of reaction channels that remove flux from the 
elastic channel at distances well beyond the strong absorp- 
tion radius, in contrast to the picture suggested by the strong 
absorption model. Continuum-Discretized Coupled-Channels 
(CDCC) calculations using either a simplified di-neutron model 

Mi 

1 ij] or a more realistic three-body model for the ^He nu- 
cleus 012-181 indicate that this long range absorption phe- 
nomenon can be explained in terms of the strong couplings to 
the breakup channels due mainly to the dipole Coulomb inter- 
action. 

Despite their success to describe the elastic data angular dis- 
tributions, it was shown in ifioll that these CDCC calculations 
failed to reproduce the energy and angular distributions of the 
a particles observed at large angles in the same experiment. On 
the other hand, it was shown in the same work that these alpha 
particles distributions could be well accounted for by assum- 
ing a two-neutron transfer mechanism leading to very excited 
states of the target. These calculations were performed using 
the distorted-wave Born approximation (DWBA), assuming a 
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two-neutron mechanism leading to both bound and unbound 
final states of the target. The relevance of the one- and two- 
neutron transfer channels has been evidenced in other reactions 
involving ^He |,2i-(5i, ,20i1 ■ 

The main diff'erence between the CDCC and DWBA calcu- 
lations is due to the couplings included in each method. This is 
schematically illustrated in Fig.[T] The CDCC method assumes 
a direct breakup mechanism in which the breakup process is 
treated as an inelastic excitation of the projectile to its con- 
tinuum spectrum (top panel). On the other hand, the DWBA 
method (or more elaborated forms of the transfer amplitude 
1 2111 ') is based on a transfer to the continuum picture, which 
treats the ^He breakup assuming a neutron transfer mechanism 
populating ^K'pb* states, as shown in the bottom panel of Fig. [T] 
Since both the projectile and target states form a complete set, 
we expect both methods to give similar results provided that 
the underlying interactions are the same and the model space is 
sufficiently large to achieve convergence of the studied observ- 
ables Ii21il . In this respect, the direct breakup and transfer to the 
continuum methods can be regarded as alternative methods to 
describe the projectile dissociation. Therefore, the failure of the 
CDCC approach to describe the a yield could be attributed to 
the truncation of the model space required in practical calcula- 
tions, and/or to the neglect of the transfer to the bound states of 
the target which, at these energies, may contribute significantly 
to the a cross section. 

The transfer to the continuum calculations presented in lfl9ll 
are based on the DWBA approximation and describe the elastic 
channel by means of a phenomenological OM, with parameters 
adjusted to the elastic data. The features of this OM potential 
must reflect the effect of channel couplings on the elastic chan- 
nel. Since the two-neutron transfer channels were found to be 
very important to describe the a yields, it would be desirable to 
investigate to what extent the explicit inclusion of these transfer 
couplings can explain simultaneously the observed features of 
the elastic scattering and the a channel, without the need of a 
phenomenological OMP adjusted to the data. To study the ef- 
fect of the transfer channels on the elastic channel one needs to 
go beyond the Born approximation, that is, to use the Coupled- 
Reaction-Channels (CRC) method. 

With this purpose in this work we present CRC and CDCC 
calculations for the ^HeH-^'^**Pb reaction at near-barrier energies, 
comparing the calculated elastic and a cross sections with the 
available data. In the CRC calculations we rely on the scheme 
illustrated in the bottom panel of Fig.IT] but we allow back cou- 
plings from the transfer to the elastic channels. We show that, 
unUke the CDCC, the CRC formalism is able to explain simul- 
taneously the long-range absorption eff'ect on the elastic scat- 
tering and the large yield of a particles for this reaction. 

Although the calculations are restricted to the ^He-i-^*'^Pb 
case, we believe that the conclusions can be extrapolated to 
other reactions induced by weakly bound projectiles. 

2. CDCC calculations 

Within the Continuum-Discretized Coupled-Channels 
(CDCC) method, the dissociation of the projectile is treated 



assuming a direct breakup picture, in which the projectile is 
excited to its unbound states (see Fig. [1). For ^He reactions, 
this has been done using a three-body reaction model (based on 
a simple di-neutron model of the projectile, '*HeH-2n ifioi [Till ) 
or a four-body reaction model (using a realistic three-body 
model of *'He='*HeH-n-Hn lfl2ilT3l[l6mil). In this work, we rely 
on the first model because this will permit a more meaningful 
comparison with the CRC calculations presented in Section [3] 
which are also based on a two-body model of the ^He nucleus. 
In particular, we will make use of the improved di-neutron 
model proposed in 111, [lot], in which the '^HeH-2n relative 
wavefunctions are calculated in a Woods-Saxon potential with 
radius R - 1 .9 fm and diffuseness parameter a^) - 0.39 fm. The 
ground state wavefunction is calculated assuming a 2S config- 
uration and an efifective separation energy of the two-neutron 
cluster Sin - 1 -6 Me V. For the a;H-^'"^Pb interaction we took 
the potential of Barnett and Lilley lf22ll . whereas the 2n-i-^*'**Pb 
interaction was approximated by the deuteron-^'^*'Pb global 
potential of Ref. [ 23.1 . ^He continuum states with relative 
angular momentum /, — 0,1,2 for the '^He-i-2n relative motion 
were considered. For the U - 2 continuum, the potential 
depth was adjusted to reproduce the known 2^ resonance 
at Ejc - 1.8 MeV above the g.s. For each value of Z,, the 
continuum was discretized using the standard binning method. 
The maximum excitation energy depended somewhat on the 
incident energy, ranging from 5 MeV (for £'iab=14 MeV) to 8 
MeV (for E\.ja-21 MeV). These calculations were performed 
with the code fresco |24]. 

The calculated elastic differential cross sections are com- 
pared in Fig. |2] (dashed lines) with the experimental data from 
Refs. |i7|-(9i]. Despite the simplified structure model used for the 
*He nucleus and the absence of any free adjustable parameter, 
the overall agreement with the data is very good at the five con- 
sidered energies. Below the barrier (£'iab=14, 16, and 18 MeV), 
the inclusion of the breakup couplings produces a decrease of 
the cross section with respect to Rutherford at cm. angles be- 
yond 60°. Above the barrier (£'iab=22 and 27 MeV) the main 
effect is the absence of the Fresnel peak. 

In Figs. [3] and 21 the calculated energy and angular distri- 
butions of the a particles (dotted lines) are compared with the 
data from Ref. 11911 (solid circles). Clearly, the calculations fail 
to describe both observables and, therefore, we conclude that 
the direct breakup model, at least within this restricted model 
space, cannot explain the large yield of a particles produced 
at large angles. On the other hand, it was shown in Ref. lfl9ll 
that these distributions could be well reproduced by means of 
DWBA calculations, assuming a two-neutron transfer mecha- 
nism. 

3. CRC calculations 

In the distorted-wave Born approximation (DWBA) calcula- 
tions of Ref. lfl9ll . the entrance channel was described using a 
phenomenological OM potential, with parameters adjusted to 
reproduce the experimental elastic angular distribution. In the 
Coupled-Reaction-Channels (CRC) calculations, we aim at ex- 
plaining the features of the elastic scattering as a consequence 
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Figure 1: Schematic representation of the couplings included in the direct 
breakup (CDCC) (top) and transfer to the continuum (DWBA, CRC) (bottom) 
calculations. 



of the same couplings which are responsible for the production 
of the a particles with the objective of reproducing both chan- 
nels within the same framework. Therefore, in order to avoid 
double counting, in the CRC calculations one has to use instead 
a bare interaction, that is, a potential that represents the inter- 
action between the colliding nuclei in absence of couplings that 
will be explicitly included. 

In particular, we use for the bare interaction the Sao Paulo 
potential (SPP) 1I25I I26I1 . This is a microscopic potential ob- 
tained by means of a double-folding procedure, using the mat- 
ter densities of the colliding nuclei and an effective nucleon- 
nucleon interaction, multiplied by an energy-dependent term, 
which accounts for part of the non-locality of the optical po- 
tential. This provides a parameter-free prescription, which has 
been successfully applied to a large variety of systems 12514321. 
Indeed, the choice of the bare interaction is not unique. The 
present choice is convenient because, being based on the mat- 
ter densities of the colliding nuclei, reduces many ambiguities 
associated with more phenomenological potentials. The bare 
interaction includes also a short-range imaginary potential of 
Woods-Saxon shape and parameters W = 50 MeV, Rq - 1 fm, 
flo^O.l fm, to simulate the in-going boundary condition for fu- 
sion. This bare interaction was recently used in the OM analysis 
performed in for the same elastic data. As it was shown in 
that work, the bare interaction alone is unable to reproduce the 
elastic data. This is illustrated in Fig.|2] where the one-channel 
calculation with the bare interaction (dotted line) is compared 
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Figure 2: Elastic scattering angular distribution as a function of the cm. scatter- 
ing angle for the reaction ''He+^o^pb at £iab = 14, 16, 18, 22, and 27 MeV. The 
dotted line is the one-channel calculation with the bare interaction, Uspp(r), the 
dashed line is the three-body CDCC calculation and the solid line is the full 
CRC calculation. The circles are the experimental data from Ref. [7|,(9(]. 



with the experimental data. At E\^i, - 22 and 27 MeV (above 
the nominal Coulomb barrier) this calculation predicts a pro- 
nounced rainbow, which is absent from the data. Below the bar- 
rier (iiiab = 14, 16, and 18 MeV), the calculation overestimates 
the data for angles beyond 50°. Thus, this prescription fails 
to describe the elastic data both below and above the Coulomb 
barrier The failure is clearly due to the influence of channel 
couplings, mainly, the 2n removal channels, which are expected 
to be very important for a loosely bound nucleus like ''He. In 
Ref. ll33ll . the effect of these channels on the elastic channel 
was taken into account adding a phenomenological component 
to the bare interaction of Woods-Saxon form. It was found that 
the real and imaginary parts of the Woods-Saxon potentials re- 
quired a large diffuseness parameter in order to reproduce the 
data. 

In the present work, we start from the same bare interac- 
tion but, instead of adding any phenomenological component, 
we introduce the two-neutron transfer couplings explicitly, us- 
ing the CRC formalism. Using the prior representation, the 
transfer couplings involve a matrix element of the operator 
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V[2n+Pb] + U[a+Pb] - ^/bare, whcre y[2n+pb] IS a real interaction 
describing two neutron states of the 2n+^"^Pb system, U[a+pb] 
is the optical potential describing the a'+^"**Pb elastic scattering 
and f/bare is the bare potential defined above, given by the Sao 
Paulo prescription. In the CRC calculations we took the same 
a+^^^Pb and 2n+^'"*Pb interactions used in the CDCC calcula- 
tions of the preceding subsection, in order to have a meaningful 
comparison between both formalisms. Nevertheless, for the lat- 
ter only the real part of the interaction was considered, in order 
to allow the inclusion of two-neutron bound states in ^'"Pb. For 
the ^He nucleus, we adopt the same di-neutron model used in 
the CDCC calculations. For the ^lOpb* final states, we consid- 
ered several values for the angular momentum of the relative 
angular motion for ^°^Pb-i-2n (If). For each value of //, the en- 
ergy spectrum was described by a set of representative states, 
following the procedure described in 1.19.1 . These representative 
states include both bound and unbound states, with respect to 
the two-neutron breakup threshold. Above this threshold, the 
2jj.208p^ continuum was discretized using 1 MeV bins, up to 
a maximum relative energy of 5 MeV. For energies below the 
threshold, we considered 5 bound states, evenly spaced in steps 
of 2 MeV. The number of partial waves If was progressively 
increased until convergence of the observables was found, in 
the kinematic region where comparison with data is possible 
(see discussion below). The coupled equations were solved it- 
erative ly using the code fresco ||24|1 . until the desired degree 
of convergence of the elastic and transfer cross sections was 
achieved. We include also the non-orthogonality correction 

The elastic differential cross sections obtained with these 
CRC calculations are displayed in Fig.|2]with solid lines. For all 
the considered energies, the agreement with the CDCC calcu- 
lation and hence with the data is significantly improved with 
respect to the single-channel calculation performed with the 
bare interaction alone. Therefore, the inclusion of the trans- 
fer channels produces essentially the same eff'ect on the elastic 
cross section as the inclusion of the ^He continuum states in the 
CDCC calculations. This is an important result which supports 
the conclusion that both, CDCC and CRC methods, populate 
to a large extent the same final states, although these states are 
expressed in different basis representations. 

Despite the good agreement between both methods, the rate 
of convergence is very different. The CRC calculations were 
performed with // - 0-8 partial waves for the 2n--'"^Pb rel- 
ative motion, whereas the three-body CDCC calculations re- 
quired only Z, = 0,1,2 partial waves for 2n-'^He in order to 
achieve convergence of the elastic scattering ifTTIl . This indi- 
cates that the elastic scattering is mostly affected by the cou- 
pling to continuum states with small energy and angular mo- 
mentum between 2n and a, and hence a representation based 
on the continuum of the projectile is more efficient to describe 
this observable. 

The calculated energy and angular distributions of the a par- 
ticles are compared with the data in Figs. |3] and |4l resp ectively. 
For comparison, the DWBA calculations of Ref. [19] are also 
included (dashed lines). We see that the CRC calculations, rep- 
resented by solid lines, reproduce the data with a similar quality 
of the DWBA calculations, although the magnitude is some- 




10 15 20 25 30 



E„ (MeV) 



Figure 3: Energy distribution of the a particles produced in the reaction 
6fjg^208p|^ at £'iab=14, 16, 18, and 22 MeV, integrated in the angular interval 
eiab = 132° - 164°. The dotted and solid lines are the CDCC and CRC calcu- 
lations. The dashed lines and circles are, respectively, the DWBA calculations 
and the experimental data from Ref. II^ . 



what underestimated. This discrepancy might be due to the 
limitations of our di-neutron model (for both the projectile and 
target states), to the choice of the underlying interactions, or to 
the contribution of other channels. Considering the simplicity 
of our model, we can say that the overall agreement is fairly 
good. Therefore, unlike DWBA calculations, the CRC calcula- 
tions are able of reproducing simultaneously the elastic and a 
production data without requiring any phenomenological fit of 
the elastic cross section. 

For a meaningful comparison with the data, we have verified 
that the number of partial waves used for the 2n-^°**Pb motion 
(// - 0-8) provides also convergence of the a cross section 
in the angular range where data exist. This is shown in Fig. |5] 
where we plot the contribution of each / / to the angle-integrated 
transfer cross section for £'iab=22 MeV. The dark histograms 
correspond to the angular range 0iab = 132° - 164°, whereas the 
dashed bars are for the full angular range. It is clear from this 
figure that the model space used in the CRC calculations yields 
convergence within the angular range of interest and hence we 
conclude that this model space is suitable to describe the a par- 
ticles at backward angles. At forward angles, the histograms 
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Figure 4: Angular distribution of tlie a particles produced in tlie reaction 
*'He+208pb at £'iab=14, 16, 18, and 22 MeV. The dotted and solid lines are the 
CDCC and CRC calculations. The dashed lines and circles are, respectively, 
the DWBA calculations and the experimental data from Ref . 1 191 . 



of Fig. |5] suggest a non-negligible contributions from // > 8 
waves. Although there are no data for the a particles emitted at 
these forward angles, it is expected that these are more easily 
described by the CDCC calculations. 

4. Trivially-equivalent local polarization potentials 

The analysis of the preceding section demonstrates that the 
two-neutron stripping channel can explain the long-range ef- 
fect found in previous phenomenological analyses of the elastic 
data L2-5]- To corroborate this conclusion in a more quantita- 
tive way, we have evaluated the so called trivially-equivalent 
local polarization potential (TELP). This is a local and L- 
independent potential constructed from the solution of the cou- 
pled equations in such a way that it simulates the net effect of 
the couplings on the elastic scattering. Here, the TELP is cal- 
culated following the prescription proposed by Thompson et 
al. ll3411 . which involves two steps. First, for each total an- 
gular momentum, a local polarization potential is calculated 
from the source term of the elastic channel equation. Then, 
a L-independent potential is constructed by averaging these 
L-dependent polarization potentials, weighted with the trans- 
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Figure 5: Two-neutron transfer cross section as a function of the 2n--'"'Pb or- 
bital angular momentum If, for the CRC calculations peiformed at an inci- 
dent energy of 22 MeV. The dashed bars and solid bars coiTespond, respec- 
tively, to the full angular range and the integrated angular distribution for 
eiab = 132° - 164°. 



fer/breakup cross section for each angular momentum. The 
TELP obtained by this procedure should be regarded as a sim- 
plified local approximation of a complicated coupled-channels 
system. If the TELP extracted from the coupled-channels cal- 
culation is a good representation of the overall effect of the cou- 
plings, the solution of the single-channel Schrodinger equation 
with the effective potential [/eft = t^bare + t^xELP should repro- 
duce the elastic scattering similar to the one obtained with the 
full coupled-channels calculation. 

In Fig. |6] we show the effective potentials (f/eif) extracted 
from the CDCC and CRC calculations at Ei^b = 22 MeV. The 
top and bottom panels correspond to the real and imaginary 
parts, respectively. The arrows indicate the radius of sensitivity 
of the real and imaginary parts of the optical potential, accord- 
ing to the OM analysis performed in |9|] . 

The dot-dashed and solid lines correspond, respectively, 
to the effective potential f/eft extracted from the CDCC and 
CRC calculations. For comparison, we have included also the 
double-folding SPP potential (dotted line in top panel), which 
corresponds to the bare interaction in the CRC calculations, 
and the phenomenological optical potential extracted in Ref. ^ 
from the fit of the elastic data (dashed line). 

For distances r > 15 fm, the real part of the effective poten- 
tial (in both the CRC and CDCC calculations) is dominated by 
a long-range attractive tail. This attractive component is known 
to arise from the dipole Coulomb interaction 133 . 35. 3( 



Al- 
though the direct breakup model provides a more natural repre- 
sentation for this effect, it is noticeable that the TELP extracted 
from the CRC shows up also this behaviour. 

At distances around the strong absorption radius the TELP 
becomes repulsive, making the effective potential shallower 
than the bare interaction. This repulsive component is mainly 
due to nuclear couplings l'33'l. The OM potential is also less in- 
tense than the bare interaction, but it does not exhibit the long- 
range tail observed in the TELPs. This might be due to the fact 
that the real part of the optical potential is mostly sensitive to 
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Figure 6: EfTective potential, f/cff=f/baic + C/telp, for ^He+^o^Pb at 22 MeV, 
extracted from the CDCC (dot-daslied lines) and CRC calculations (solid lines). 
The dashed line is the phenomenological optical model potential extracted in 
Ref. 1 9] from the fit of the elastic data. The dotted line is the bare interaction 
used in the CRC calculations. The top and bottom panels correspond to the 
real and imaginary parts, respectively. The inset in the top figure compares the 
polarization potential 1/telp extracted from the CDCC and CRC calculations 
with the bare interaction. 



distances around the strong absorption radius. 

In the bottom panel of Fig. |6] we see that the imaginary 
part of the effective potential extracted from both, the CDCC 
and CRC calculations, is absorptive and exhibits also a diffuse 
tail. This behaviour is also observed in the optical potential 
(dashed line) and is related to the long-range absorption effect 
discussed in previous optical model analyses of this reaction 

13- 

1^ . Therefore, another important conclusion of this work is that 
the CDCC and CRC calculations provide a microscopic inter- 
pretation of this long-range absorption effect in terms of chan- 
nel couplings. Quantitatively, the effective potential extracted 
from the CDCC calculations is closer to the phenomenological 
optical potential. We emphasize, however, that due to the L av- 
erage, the TELP is not strictly equivalent to the full set of equa- 
tions, and hence only the qualitative behavior is meaningful. In 
addition, the transfer couplings are intrinsically non-local and 
thus their approximation by a local object has to be interpreted 
with caution. 



5. Summary and conclusions 

We have studied the elastic and a-production channel for 
the ^He-i-^'"^Pb reaction at energies around the Coulomb bar- 
rier. The available experimental data have been compared with 
three-body Continuum-Discretized Coupled-Channels (CDCC) 
and Coupled-Reaction-Channels (CRC) calculations using in 
both cases a di-neutron model of the ^He nucleus. 

In the CDCC calculations, the projectile dissociation is taken 
into account by including the coupling with the "*He-i-2n con- 
tinuum states. The inclusion of these couplings have a signifi- 
cant effect on the calculated elastic cross section, and provides 
a good description of the experimental elastic data. However, 
the calculated a cross sections largely underestimate the exper- 
imental data, in agreement with the findings of Ref. 111. 

In the CRC calculations, we consider the coupling to 2n 
transfer channels and its influence on the elastic scattering cross 
section. For the bare interaction, we have used the microscopic 
double-folding Sao Paulo potential, supplemented with a short- 
range Woods-Saxon imaginary potential to account for com- 
plete fusion. The inclusion of two-neutron transfer channels 
produces a strong effect on the elastic cross section, providing 
a very good agreement with the experimental data. Moreover, 
these calculations explain also reasonably well the energy and 
angular distribution of the a particles measured at backward an- 
gles. In particular, the shape of the energy distribution is very 
well accounted for, although the magnitude of the cross section 
is somewhat underestimated. This discrepancy might be due to 
the limitations of our simple di-neutron model or the choice of 
the interactions. In addition, besides the two-neutron transfer 
channels, other channels could contribute to this reaction. For 
example, in Ref. i37tl the authors studied the effect of the one- 
neutron stripping channel, (^He,^He), on the elastic and fusion 
cross section for several targets. In the ^He-i-^°^Pb case, they 
found that the inclusion of these channels produce a reduction 
of the elastic cross section in the region of the Coulomb rain- 
bow, and an increase at backward angles. The simultaneous 
inclusion of both mechanisms would be of interest. 

The trivially-equivalent local polarization potentials (TELP) 
derived from the CRC and CDCC calculations exhibit the char- 
acteristic long-range real and absorptive parts, in agreement 
with the behavior observed in the phenomenological optical po- 
tentials extracted in previous OM analyses of the same data. 
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